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The coupling of distinct systems underlies nearly all physical phenomena and their applications. A 
basic instance is that of interacting harmonic oscillators, which gives rise to, for example, the phonon 
eigenmodes in a crystal lattice. Particularly important are the interactions in hybrid quantum 
systems consisting of different kinds of degrees of freedom. These assemblies can combine the 
benefits of each in future quantum technologies. Here, we investigate a hybrid optomechanical 
system having three degrees of freedom, consisting of a microwave cavity and two micromechanical 
beams with closely spaced frequencies around 32 MHz and no direct interaction. We record the first 
evidence of tripartite optomechanical mixing, implying that the eigenmodes are combinations of one 
photonic and two phononic modes. We identify an asymmetric dark mode having a long lifetime. 
Simultaneously, we operate the nearly macroscopic mechanical modes close to the motional quantum 
ground state, down to 1.8 thermal quanta, achieved by back-action cooling. These results constitute 
an important advance towards engineering entangled motional states. 



One of the present goals in physics is the explanation of 
macroscopic phenomena as emerging from the quantum- 
mechanical laws governing nature on a microscopic scale. 
This understanding is also important for future quantum 
information applications [iHj], since many of the most 
promising platforms base on nearly macroscopic degrees 
of freedom. For macroscopic mechanical objects, their 
potential quantum behavior [f| has been actively inves- 
tigated with resonators interacting with an electromag- 
netic cavity mode @. In particular, freezing of the me- 
chanical Brownian motion 0, HJ down to below one quan- 
tum of energy has recently been observed 0, . 

Another result of radiation-pressure interaction is the 
mixing of the normal modes of vibration into linear com- 
binations of the uncoupled phonon and photon modes. 
This has been verified with optomechanical Fabry-Perot 
cavity 11 1 , and by its on-chip microwave analogy [t| us- 
ing an aluminum membrane. All this work is paving the 
way towards engineering non-classical motional [H[ and 
hybrid quantum states [l], [!, 0, [HJ for basic tests of quan- 
tum theory [l3 , 14 1 , as well as applications in foreseeable 
future. 

Corresponding phenomena in systems comprising of 
more than two active degrees of freedom [TBI, [l6| have 
received less attention. In the optomechanical crystal 
setup, coupling of mechanical vibrations via radiation 
pressure interaction has been demonstrated with the zip- 
per cavity [17| . however, with some direct interaction be- 
tween the beams. In the microwave regime, such mea- 
surements have remained outstanding. In this work, we 
take a step further, and examine a multimode system 
where two micromechanical beams, having resonant fre- 
quencies wi/(2tt) = 32.1 MHz and w 2 /(27r) = 32.5 MHz, 
are each coupled to a microwave on-chip cavity. We ob- 
tain the first evidence of hybridization of all the three de- 
grees of freedom. This was made possible by operating, 
as opposed to the optical setup, in the good-cavity limit, 
that is, both wi and L02 were much larger than the cavity 



decay rate 7 C . Simultaneously, the mechanical modes are 
found at low occupation numbers n m near the quantum 
ground state, down to n m = 1.8 ± 0.5. These are the 
lowest occupation recorded with nanowire resonators to 
date, and they occur via the sideband cooling, starting 
from operation at dilution refrigerator temperatures at a 
few tens of mK. 



In the optical domain [17j . a three-mode optome 



chanical system can be pictured as a Fabry-Perot cav- 
ity with both massive mirrors mounted via springs, see 
Fig. la. In the microwave version, (Fig. lb), the mir- 
rors are replaced by movable capacitors formed by the 
mechanical resonators. The electromechanical coupling 
arises when each mechanical resonator (labelled with j), 
expressed as time- varying capacitors Cj, independently 
modulates the total capacitance, and hence, the cavity 
frequency oj c . This is described by the coupling energies 
g.j = (uj c /2C)dCj/dxj. 

In the actual device (see micrograph in Fig. lc), the 
mechanical resonators are spatially separated by about 
100 microns, and we hence expect the direct interaction 
between the beams via the substrate material to be negli- 
gible. We drive the cavity strongly by a microwave pump 
signal applied at the frequency up near the cavity fre- 
quency. The pump induces a large cavity field with the 
number of pump photons np ^> 1. This allows for a lin- 
earized description of the electromechanical interaction 
and results in substantially enhanced effective couplings 



Results 

Cavity with several mechanical modes The gen- 
eral Hamiltonian for ./V mechanical resonators with fre- 
quencies ujj, individually coupled to a common cavity 
mode, written in a frame rotating with the pump, be- 
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FIG. 1: Cavity-controlled coupling of micromechanical 
resonators, a, Optomechanical system with both mirrors 
movable has three degrees of freedom, b, analogy of (a) with 
a microwave-regime cavity represented by the lumped element 
inductors and capacitors and arbitrary number j = 1...N of 
mechanical resonators. The experiments are performed by 
shining a strong pump microwave tone (frequency up), and 
possibly a weak probe tone (frequency Ui n ) in a dilution re- 
frigerator down to 25 mK temperature, c, images of the all- 
aluminum superconducting sample, showing four fabrication 
jigs for a beam. Total of three beams (X1...X3) were oper- 
ating. Beams 1 and 2 had large electromechanical couplings 
gi/2% = 1.8 MHz/nm, g2/2n = 2.0 MHz/nm, whereas the 
third beam had an order of magnitude smaller coupling. 



comes 

JV N 

H = -ftAata + h 0Jjb]bj - h (at + a) ^ Gj [b] + b 

(1) 

Here, a is the annihilation operator for the cavity mode, 
bj are those of the mechanical resonators, and xoj = 
jUj is the zero-point amplitude. If one assumes 
roughly similar mechanical resonators, viz. u)j ~ u) m 
and Gj ~ G, the cavity becomes nearly resonant to 
all of them at the red-sideband detuned pump condition 
A = -uj m . 

Coupling of the mechanical resonators via the cavity 
"bus" can be anticipated to be significant if the me- 
chanical spectra begin to overlap when their width in- 



creased by the radiation pressure, 7 e g = 7 m +7opti grows 
comparable to the mechanical frequency spacing. Here, 
7opt = 4G 2 /7 C . The equations of motion following from 
Eq. (QJ allow one to verify the above assumption. We 
will in the following focus on two resonators, N = 2. 
The result of such calculation, which is detailed in the 
Supplementary Information, is that they experience an 
effective coupling with energy Gi 2 = G\Gihc- 

In the limit of strong coupling, \J(G\ + G\)j2 > 7 c /4, 
the system is best described as a combination of two 
"bright" modes with linewidths ~ 7c/2, and of a "dark" 
mode having a long lifetime ~ 7 m . In the latter, the two 
mechanical resonators oscillate out-of-phase and synchro- 
nize into a common frequency (cji + 012) /2. In any case, 
the mode structure can be written with the normal co- 
ordinates for mode n, for position quadrature according 



to = D ( c n) q c 



and similarly for the 



momentum quadrature P( r 

In the experiment, we use doubly clamped flexural 
beams which are made by the use of ion beam milling 
of aluminum (lil ]. having an ultranarrow ~ 10 nm vac- 
uum slit to the opposite end of the cavity for maximizing 
the coupling. The cavity is A/2 floating microstrip sim- 
ilar to Ref. resonating at lo c /(2-k) = 6.98 GHz. 
The total cavity linewidth 7 C = 7^ + 7/ ~ (27r) x 6.2 
MHz is a sum of the internal damping 7//(27r) = 1.4 
MHz, and the external damping ^ e /(2tt) = 4.8 MHz. At 
at the highest powers discussed, however, we obtain a 
decreased j]/(2ir) — 0.9 MHz, typical of dielectric loss 
mechanism. 

There are a total of three beams as shown in Fig. lc. 
Two of the beams have a large zero-point coupling of 
giXoi/2-K = 39 Hz, and 32^02/2^ = 44 Hz. The frequen- 
cies of beams 1 and 2 were relatively close to each other, 
u) 2 — ijJ\ ~ (2ir) 450 kHz, such that it is straightforward to 
obtain an effective coupling G12 of the same order. The 
third beam had an order of magnitude smaller coupling, 
and hence we neglect it in the full calculation, setting 
N = 2 in Eq. (JTJ in what follows. 

The measurements are conducted in a cryogenic tem- 
perature below the superconducting transition tempera- 
ture of the Al film (~ 1 K) in order to reduce losses. We 
used a dilution refrigerator setup as in Ref. [3] down to 
25 mK, which allows for a low mechanical mode temper- 
ature in equilibrium. Apart from this, cryogenic temper- 
atures or superconductivity are not necessary (20j . 

Experimental data The modes can be probed by 
adding to the experimental setup another, the probe, 



tone 



11, m 



2l| with frequency Ui n . The maximum 
pump power we can reach is 4 /iW, inducing a coher- 
ent np ~ 6 x 10 8 photon occupancy in the cavity, and 
effective coupling G 12 = (2ir) x 150 kHz (with Gi = 0.93 
MHz, G 2 = 1.0 MHz). In Fig. 2a, we display changes 
of the cavity absorption by increasing the pump power. 
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FIG. 2: Hybridization of microwave photons and two radio-frequency phonons. a, Measured reflection coefficient of 
probe microwave while a strong pump is applied at the red sideband, A ~ —cji, while increasing pump amplitude, from bottom 
to top. b, c, zoomed-in view of (a). An excellent agreement with the theory is obtained with reflection calculation using a 
complete model based on Eq. ([1]), for two beams 1 and 2 (solid lines). The dashed lines show a fit with an independent-resonator 
model for beam 1. The narrow feature at 6.982 GHz is due to beam 3. The curves are displaced vertically by 3 dB for clarity, d, 
normal modes of oscillations expressed as linear combinations of photon and two phonons for the dark mode. The parameters 
are as in c. The arrow denotes the highest experimental value, e, Theory plot for higher up; the dark mode is expected to 
become narrower and more and more phonon-like. Here, A = — (wi + W2)/2. For simplicity, we take g^ = in d and e. 



The two peaks at the bottom of the cavity dip are due 
to the beams 1 and 2. They broaden with increasing np, 
finally leaving behind a narrow dip between them. The 
overall lineshape is clearly not a sum of two Lorentzian 
peaks. 

In Fig. 2c we show a zoom-in view of the two peaks, 
together with theoretical predictions. One immediately 
sees that an attempt to simulate either peak by a single 
mechanical resonator coupled to cavity (N = 1 in the 
analysis), thereby neglecting the cavity- mediated cou- 
pling, fails to explain the lineshape at large np > 10 8 . 
The full simulation with both beams 1 and 2 included, 
however, produces a remarkable agreement to the exper- 
iment. In order to create these curves, we further take 
into account a shift of the mechanical frequencies due to 
second-order interaction given as —\g^ n P, where 

gf = x 2 0j u c /2C (d 2 C ] /dx 2 j ) ~ (2tt) 3 x 1(T 4 Hz is the 
second-order coupling coefficient. The third beam, visi- 
ble in Fig. 2c as a narrow feature just below 6.982 GHz, 
was later fitted by running a single-resonator calculation. 
This is justified because its mixing is negligible owing to 
weak coupling. 

The narrow dip between the peaks manifests the on- 
set of the dark mode, where the cavity participation ap- 
proaches zero with growing Gj , as occurs also in other tri- 
partite systems [22|. As usual for weakly decaying states, 
the dark modes can be useful for storage of quantum in- 
formation 17|, l23[. The other two dips (best discerned 



in the theoretical plots at higher Gj, sec Fig. 2c) are the 
bright modes, and they retain fully mixed character of all 
three degrees of freedom. In Fig. 2d we display the pre- 
diction for the mode expansion coefficients with respect 
to coupling. Asymmetry in the role of the mechanical 
modes is sensitive to parameters. With the maximum 
coupling in the experiment (Gj/j c ~ 0.15), the normal 
modes are well mixed: the dark mode can be expressed 
as [D C ,D 1 ,D 2 ] ~ [0.3,0.3,0.9]. 
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FIG. 3: Pump detuning measurements, a, Measured 
spectrum of the tripartite cavity electromechanical system 
when pumped at a fixed input power such that at resonance, 
A ~ (— u>i, — 0J2), np ~ 5.6 x 10 8 . The pump frequency was 
slowly swept from left to right, b, theory prediction. 



One can also vary the pump detuning, thereby effee- 
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FIG. 4: Sideband cooling of motion down to 1.8 thermal quanta, a, cavity output spectra representing the thermal 
motion peaks of beams 1 and 2 measured at 345 mK and 35 mK, with np = 4 x 10 J . b, calibration of the detection according to 
temperature dependence of the mechanical energy n m hhJ m = ksT, given by the area of the peak (data for beam 2). c, output 
spectra of beam 2 at pump occupation np increasing from top to bottom. The curves are shifted vertically by 0.3 y-axis units 
for clarity. The corresponding mechanical occupations are written on the right, d, progress of back-action cooling (beam 2) as 
a function of pump occupation. The mechanical mode temperature is given by circles, and the cavity temperature by squares. 
Red: high temperature, 245 mK; blue, base temperature (28 mK). The ideal theoretical behavior are plotted by the solid lines. 
Unless the data point for njT, including its error bars, was of positive sign, only the absolute value of the error bar is plotted. 



tively detuning the cavity and the mechanical resonators 
111 . l2l| . An anticrossing of the cavity and mechanical 
modes is seen in Fig. 3a up to detuning A ~ — oj m , above 
which the system exhibits parametric oscillations. The 
pertinent simulation, Fig. 3b, portrays the main features 
of the measurement, except bending of the cavity fre- 
quency towards lower value in an upward sweep of A. 
We expect such possibly hysteretic behavior of the cav- 
ity frequency to be due to second-order effects beyond 
the present linear model. 

Sideband cooling We now turn the discussion into 
showing how the tripartite system resides nearly in a pure 
quantum state during the mode mixing experiments. 
The mechanical resonators can be sideband-cooled with 
the radiation pressure back-action 



24-261 under the 



effect of the pump which also creates the mode mix- 
ing. Although the ground state of one mode was re- 
cently achieved this way 0, , cooling in multimode or 
nanowire systems is less advanced. 

The incoherent output spectrum can be due to cither 
a finite phonon number n m , or non-equilibrium popula- 



tion number n<T of the cavity. The former manifests as a 
Lorentzian centered at the upper sideband, as in Fig. 4a, 
with the linewidth y e g = j m + AG 2 /j c , whereas the latter 
gives a small emission at broader bandwidth. The ana- 
lytical theory [H and experimental details are discussed 
in Supplementary Information. 

At low input power, when the cavity back-action 
damping is small and ~ 0, phonon number is set 
by thermal equipartition n m jhujj = ksT. This is best 
observed by the area of the mechanical peak in the out- 
put spectrum, as in Fig. 4a, which should be linear in 
cryostat temperature. The linearity holds down to about 
150 mK temperatures (Fig. 4b), below which we observe 
intermittent heating which is sensitive to parameters. 

The resulting cooling of beam 2 is displayed in Fig. 4d. 
At a relatively high temperature of 245 mK, the data 
points follow theory well up to input powers correspond- 
ing to np ~ 5 x 10 7 . At the minimum cryostat tempera- 
ture of 28 mK, the mode is not thermalized in a manner 
which depends irregularly on power. We obtain a mini- 
mum phonon number n m i ~ 1.8, where the mechanical 
mode spends one third of the time in the quantum ground 



5 



state. The other mechanical mode, beam 1, cooled simul- 
taneously down to 2.5 quanta, possibly compromised by 
the slightly smaller coupling. The cooling is quite clearly 
bottlenecked by heating, by the pump microwave, of the 
bath to which the mechanical mode is coupled. For the 
optimum data points, the starting temperature for cool- 
ing has raised to 150 mK, and above this it grows roughly 
quadratically with np. 



Discussion 

The coupling of micromechanical resonators mediated 
via microwave photons in an on-chip cavity is a basic 
demonstration of the control of a multimode mechani- 
cal system near the quantum limit. The setup provides 
a flexible platform for creation and studies of nonclassi- 
cal motional states entangled over the chip j27|, or over 
macroscopic distances 28, [2!|. The setup is easily ex- 
tended to embody nearly arbitrarily many mechanical 
resonators, hence allowing for designing an electrome- 
chanical mctamaterial with microwave-tunable proper- 
ties. For future quantum technology applications even 
at elevated temperatures [20[, data may be stored [Io| 
in weakly decaying dark states identified in the present 
work. 
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Hamiltonian 



Methods 

Device fabrication 

The lithography to define the meandering cavity and jigs 
for the beams consisted of a single layer of electron-beam ex- 
posure, followed by evaporation of 150 nm of aluminum on 
a fused silica substrate. Since the cavity is of planar struc- 
ture, it has no potentially problematic cross-overs which tend 
to have weak spots in superconducting properties and hence 
could heat up the cavity. We thus expect the cavity to portray 
a high critical current limited only by the intrinsic properties 
of the strip. With a 2 (im wide strip, we expect the criti- 
cal current to be several mA. With the maximum circulating 
currents in cavity (Fig. 2), about 0.5 mA of peak value, there 
was no sign of nonlinearity in the cavity response. 

The mechanical resonators were defined by Focused Ion 
Beam (FIB) cutting, as in Ref. [HI ]. We used low gallium 
ion currents of 1.5 pA which gives the nominal beam width 
of about 7 nm. In order to maximize sputtering yield with 
minimal gallium contamination, we used a single cutting pass 
mode. Otherwise, the cut beams tend to collapse to the gate. 
With the current recipe, fabrication of down to 10 nm gap 
widths over 5-10 microns distance can be done with about 50 
% yield. 

Theory of cavity-coupled resonators 

The Hamiltonian for N mechanical resonators each coupled 
to one cavity mode via the radiation pressure interaction is 

JV JV 

H — hjj c a?a + h^^cjjb^bj — fog) a^~] ffj^Qj (b] + bjj + Hp 

(2) 



2hy P £J E (a 1 + a) coscjpt induces a large cavity field 
with the number of pump photons np 3> 1. 

Under the strong driving, the cavity-mechanics interac- 
tion can be linearized individually for each beam, resulting 
in Eq. JT|). Further details of the calculations are given in the 
Supplementary Information. 

Sideband cooling 

For one beam, the measured output spectrum divided by 
system gain, S ut{w), carries information on the phonon num- 
ber according to 



„ _ 7b T , IE 

>3out — n c i 7opt~ — 

7c 7c (w 



7off 



) 2 + 7 c V4 



(n m - 2n^ 



(3) 



3=1 



which is a Lorentzian centered at the cavity frequency. A non- 
zero base level is due to possible broadband emission from the 
cavity, due to a thermal state with occupation n^ . We sup- 
pose that at modest pump (cooling) powers of utmost interest, 
coupling of the mechanical resonators remains insignificant, 
and thus we can apply Eq. (0 separately for each resonator. 
Some deviations can, however, occur at the highest powers 
with n P ~ 3 x 10 s . 

In order to obtain the phonon number under strong back- 
action, we fit a Lorentzian to each peak as in Fig. 4c, obtaining 
independently the base level, amplitude and linewidth for ev- 
ery input power. These values are then compared to Eq. <(3j . 
This leaves, however, yet too many unknowns in order to ob- 
tain the occupations of the mechanics and cavity. This is 
basically because of uncertainties in the attenuations of both 
input and amplifier lines, which limit the accuracy of estimat- 
ing np. We can obtain a calibration using the linear temper- 
ature dependence, and by fitting the dependence of 7 c ff on 
input power. 

As opposed to for instance to Ref. the bottleneck for 
cooling is not heating of the cavity by strong pump, as seen in 
Fig. 4d, where the cavity temperature is increasing only little 
up to the optimum cooling powers corresponding to about 
0.15 mA of peak current values in cavity. This can be due to 
the simplistic structure of the cavity. 

The above analysis supposes the presence of nearly no non- 
equilibrium photons in the cavity at zero or low input powers. 
This situation can be investigated by observing a possible 
emission from about cavity linewidth under these conditions. 
A complication arises because the cavity linewidth is so large 
such that a small change is easily overwhelmed by modest 
standing waves in cabling. We avoided this issue by using 
temperature dependence of the cavity frequency, due to the 
kinetic inductance in the long microstrip. Thereby, a broad 
emission peak moving according to the cavity frequency would 
be easily distinguishable. We observed no such signal down 
to the level of rip < 0.1 which justifies the above analysis. 
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RESPONSE OF COUPLED RESONATORS 

We detail here the theory of the response of a driven cavity coupled to N (with emphasis on the case N = 2) 
mechanical resonators. In particular, we show how it is possible to describe the system from two complementary 
points of view. On one hand, the cavity-mechanics interaction provides a dressing of the dynamics of the resonators, 
leading to the definition of effective mechanical frequencies and dampings. On the other hand, for strong coupling, 
we find that on top of the normal mode splitting found for single resonators, the system exhibits a dark mode which 
gets asymptotically decoupled from the cavity, and whose linewidth tends to the bare linewidth of the mechanical 
resonances. 

We define the following symbols: The mechanical resonators have frequencies Uj and linewidths jj. Their annihila- 
tion operators are denoted by bj while qj — (b^ + bj^j and pj = ^= (b^ — b^j arc the dimcnsionlcss amplitude and 

momentum of the mechanical oscillations. The mechanical zero-point amplitude is x j = y/H/2mjUJj. Correspond- 
ingly, the cavity has the frequency uj c and linewidth j c ~ +"fi, consisting of the external and internal dissipation, 
respectively, a is the annihilation operator for the cavity. The pump microwave (applied at the frequency cop near 
the cavity frequency) has the detuning from the cavity A = up — oj c , and the power Pp. 

The Hamiltonian for N mechanical resonators each coupled to one cavity mode via the radiation pressure interaction 

is 

N N 

H = tkj c a^a + h^^uijbjbj — ha^a gjXaj (bj + bj^j + Hp (4) 
j=i j=i 



The pump with the Hamiltonian Hp = 2h^J P £J E (af + a) coscjpt induces a large cavity field with the number of 
pump photons np 3> 1 . The number of quanta the in cavity due to pump is 

n P = — 5- (5) 

where angular frequency units are used. 

The equations of motion following from the Hamiltonian, Eq. (1) in the main text, linearized about a steady-state, 



can be written as [31 [ 



N 

a = iAa - -j a + i Yl G i ( H + V^i a L ( 6 ) 

3 = 1 i=I, E 

q 3 = UjPj (7) 

Pj = - u 3<l3 - IjPj + G 3 (V + a) + (8) 

where a\ n are input fields to the cavity (/ and E denoting internal and external input fields, the previous describing 
internal dissipation and the latter the coupling to the measurement setup), and £j are fields driving the mechanical 
resonators. Moreover, Gj = \Jnpg 3 - describes the cavity driving (by photon number np) enhanced coupling between 
the cavity field and resonator j. Fourier transforming and, after some algebra, assuming a probe field ai n = af^ and 
neglecting the noise terms £j yields 

N 

X = a) +a = x(u))Y2, G j<lj - i \ / TE X-A M a in + xa M (9) 

3 = 1 

q j =Xj(«>)G j X (10) 
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where the response functions are 



and 



Solving for X, we get 



with 



Xa(w) = -r 1 . /0 , Xi( w ) = -3 ^— (11) 



2A 

X(w) =X-a(w) -XaH = 7 — — — ttt^ r^. (12) 



X = y/jE~M(uj) X-Aain + XAO-ln (13) 



M(u) = jv . (14) 

This field is coupled to the cavity output field by 

a& k = y/7EaW-a& ) . (15) 



When the effect of a\ n on a can be disregarded (see below), we can write a out = R(u))a- m , where R(uj) = 
— ijeM ( w )x~a(w) — 1 is the reflection probability amplitude in the cavity. Figure 2 of the main text shows its 
absolute value, corresponding to the measured observable. 

The effective dynamics of the system is qualitatively different in different regimes. For typical structures -f c ^> -fj. 
In this case, when G\ + G\ <C (-f c /A) 2 , the cavity response x( w ) i s approximatively independent of frequency within 
the relevant frequencies describing the mechanics. In this case, we can view the system as two mechanical resonators 
effectively coupled by the cavity (Sec. J. On the other hand, for large Gj (sec. IB), the normal modes of the total 



system split. There are two bright modes strongly coupled to the cavity, with resonance frequencies cj& ps A± yX^j @j 

and linewidths ~ 7c/2, and one "dark" mode at around (uj\ + u)%)/2 and with a linewidth ~ jj. The latter is almost 
decoupled from the cavity, with the coupling decreasing with an increasing Gj. 



Weak-coupling regime 

In the following, we concentrate on the case of red-detuned driving, A rj — uij in the sideband resolved case 7 C -C u>j, 
and study the response at frequencies u> ~ — A. This allows us to approximate 

XiH ~ ir, -x — , X~ — , a _l • — M- ( 16 ) 

2(ujj — ui) — I'fj uj + A + z7c/2 

The same approximation yields xa(^) *C x-a(w) and we may hence disregard the input field a in at these frequencies. 

Let us restrict to the case of two mechanical resonators. Then the response function M(u>) may be written in the 
form 

M(ui) = (wi - w - i7i/2) (w 2 - w - 172/2) 



(ui - uj - iyi/2) (0J2-OJ- 272/2) - xM [Gj (uj 2 -uj- 172/2) + G 2 2 (wi - uj - «7i/2)] /2 
_ (wi +u) — i-fi/2) (uj 2 +u> — 472/2) 
= [w{ + w - z 7l c /2) (uj c 2 +uj- i l2 ]/2) 

We return to the definition of uij and 7? below, but now we rewrite the response in a more compact form 



(17) 



X = — — X-AV^*in (18) 
Xl X2 

x c 

Qi = G\— X\X-A^flEa- m (19) 

X2 

Xi 

92 = G 2 X2X-A\/7B a in, (20) 

Xi 
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where Xj — Xji^j ^ This form is readily generalized for more than two mechanical resonators. We now 

calculate uj^ and 7?. For the sake of compactness, we include the broadening to the frequencies, allowing them to 
assume complex values 



i 7i /2. 



(21) 



The frequencies wj, uj 2 appearing in the denominator of Eq. (|17[) can be expressed as 



(22) 



where 



w? = Wj - x(w)Gj/2 



(23) 



are the effective frequencies (and dampings) of the mechanical resonators induced by the presence of the cavity, in 
the absence of the other mechanical resonator. The frequencies (and dampings) lo? are thus the dressed mechanical 
frequencies in the presence of the cavity-mediated coupling between the mechanical resonators. Equation (|22j) describes 
an avoided crossing of the two mechanical resonators coupled by the energy 



2 7c 



(24) 



The latter form is valid in the limit lo c ^ + A <C 7c , where we can approximate x ~ —2i/j c . For strong coupling or 
large initial frequency difference, the determination of io\ 2 requires the solution of Eq. (|22p with the replacement 
x(w) — > x( w i 2)- The solutions of such a calculation are shown in Fig. [5j plotting the resulting frequency shifts and 
changes in the damping. 

When G12 <C uj\ — wf > we can disregard the coupling between the mechanical resonators, and the response of the 
cavity is a simple product of the individual responses. In the experiments, we reach G\ 2 = (2ir) x 150 kHz (with 
G\ = 0.93 MHz, G2 = 1.0 MHz), which is roughly one third of the frequency difference between the mechanical 
resonators, and allows for the coupling effect to show up in the shape of the response curves. This is shown in Fig. [51 
which compares the coupled response to the case of individual resonators. 

Although the dressed response is affected by the coupling, the frequency shifts are not easily observed from the 
cavity response as they tend to be overwhelmed by the growing width ~ Im(w|). However, the frequency shifts can 
be seen in the mechanical response q{ui)/a- m , plotted in Fig. [7] 



Strong-coupling regime 



In the strong-coupling limit xG\G\ > (u)f 



<\2 



lo\ 2 can be expressed as 



J l,2 



oj 1 +LJ 2 -x{Gi + G 2 ) /2 ± a/x 2 {G\ + G%) /4 



from which we obtain two equations for and 



bj.-, 



2 

Ll>1 + 0J 2 



+ XK) (Gl + Gl) 12 



(25) 

(26) 
(27) 



The second equation defines an asymptotic dark mode residing at the frequency between the two mechanical frequen- 
cies, and having a damping (71 + 72V2, i.e., typically much lower than that of the other dressed modes. For the first 
equation we have to include the full frequency dependence of x{uj) from Eq. (|16[) . Multiplying the equation by x' -1 ' 
then yields another second-order equation for wj, with the roots 



uil = - ( -2A + wi + uj 2 + i"i c ± V (-i7c + 2A + uj x + u 2 ) 2 + 8Gj + 8G 



(28) 



10 




FIG. 5: Values of ui\, lo 1 ^ in the weak-coupling regime. Left: real part of the frequencies and right: effective damping. On 
the left, the dashed lines show the approximation with a frequency independent coupling Gvi and disregarding the optical 
spring effect (Eq. p3[) ), whereas on the right, the dashed line shows the imaginary part of the effective frequencies wj. For the 
damping, the coupling between the resonators shows up only in the strong-coupling regime. Here and in Figs. 17181 the curves 
have been calculated with A = 57 c , wi,2 = |A| ± 0.0075A and 71 = 72 = 3 x 10 -5 |A|, close to the experimental values, and for 
simplicity Assuming G\ — G2. 




u/|A| u/|A| 

FIG. 6: (left)Comparison of the reflection probability \R(yS) | for a single oscillator vs. that of two oscillators, (right) Comparison 
of the reflection probabilities when including (solid line, calculated with ujj) or excluding (dashed line, with ojJ) the coupling 
induced normalization of the frequencies (G = 0.l7 c ). 





8 . 1 (2 res. case) 




S. 1 (1 res. case) 




FIG. 7: Mechanical spectrum of \qi\ in the case of two (blue line) and one (red line) mechanical resonators, showing the 
frequency shift due to the coupling with the other resonator, an analogous situation would be observed for I52I (G — 0.l7 c ). 
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For simplicity, we consider in the following the case A = — (lu\ + uj 2 )/2 which simplifies the expression to 



= -A + \ U 7c ± ^8(G? + G§) - 7c 2 J ■ (29) 

When G\ + G\ > 7c /8, the real part of the frequencies tends to 

ul -A ± G, (30) 



where G = \J (G\ + G\)/2. Moreover, the linewidth of these modes, whenever G > 7 c /4, is given by 7 c /2. In this case 
these frequencies and linewidths can be seen in the cavity response, as the cavity response function |i?(w)| has three 
dips: the outer dips, of width ~ 7c/2 correspond to the (bright) modes wj, whereas the inner narrow dip corresponds 
to the (dark) mode W2 (see Fig. 2e of the main text). 

Let us try to understand these modes from starting the equations of motion by redefining the mechanical motion 
to relative and center of mass motion (weighted by the couplings Gj), i.e., 

qa = (G 2 qi-G 1 q 2 )/(G 2 + G 1 ) (31) 
q s = (G iqi + G 2 q 2 )/{G 2 +G 1 ). (32) 

Below, we show that in the strong-coupling limit the above frequency a;!; describes the "dark" mode q a and the two 
frequencies ui1± the "bright" modes q s . 

The equations of motion for these modes are 

a = iAa + i(Gi + G 2 )q s + y/jEai n (33a) 

q s = UJzPs + UJAPa (33b) 
<la = W S p a + UJAPs (33c) 

p s = -UJsq s + u} A q a - IPs + J^ lG % (a 1 + a) (33d) 

Cxi + (j 2 

Pa = -^S<7a + W A <7s - IPa (33e) 

with aj£ = (loi + uj 2 ) /2 and a; a = (^i —^2) /2. We have assumed that the linewidths of the two mechanical 
resonators are the same. If there is no frequency mismatch, the symmetric and antisymmetric modes are uncoupled. 

Solving Eqs. (f3"3"]l . or substituting Eqs. (|19l20p into (j3Tj) (and using the complex frequency notation) we obtain, in 
the strong-coupling limit 

4GiG 2 — 
qs = tr Z7M 7 rT7 TTT jTV7^ a m (34) 

4G 1 G 2 uj a . — ,„_.. 

[G 2 + Gi) [w — (w — J (^ — ^> 2 ) 

Considering that = ws, and assuming for simplicity G\ = G 2 = G, we have, when u ~ &i±, 



Vis , 

Qs+ = ( TT a m 36) 



(u -oj1 + ) 



the cavity field becomes 



Analogously, when u> ~ — A, we get for <7 a 



x = or , ,c ( 38 ) 
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FIG. 8: Comparison between the calculated values of wj^ with the values obtained by diagonalization of the equations of 
motion (solid lines). (Left:) Real part of the frequency. Dashed line is the weak-coupling expression obtained from Eq. ()22[) by 
using Eq. (|24[1 and the dash-dotted line is obtained from the strong coupling expression Eq. (|28[1 , (Right): Effective linewidths 
of the modes. Black dashed line is from eqs. (|23[) . showing the regular optically induced damping, whereas the dash-dotted 
lines are obtained from Eq. (|28l) . These frequencies have been calculated with the same parameters as in the previous figures. 



In this case the cavity field around the dark mode resonance frequency becomes 



(wi + A) (u)2 + A) — 
X = ^— ; \/jEa ir 



G 2 (A - uj\) 



(40) 



From Eq. (|39[) we see that, in the strong coupling limit, the mechanical mode q a is asymptotically decoupled from 
the input as well as from the cavity field, allowing to identify this mode as a dark mode. 

From a somewhat different perspective, the above picture and the values of the effective frequencies can be recovered 
by diagonalizing the equations of motion for q ai p a , q s , p s , and a. In FigO we have compared the relevant eigenfre- 
quencies around — A as obtained by the diagonalization of Eqs. ()33[) with the strong and weak coupling expansion for 
w 1± and u>2- In the weak-coupling regime the cigenstates correspond (in the limit — > 0) to the normal modes of 
the uncoupled system (i.e. cavity, mechanical resonator 1, mechanical resonator 2, cavity field). The corresponding 
cigenfrequencies match the expression given in Eq. (|22|) in the case where x is independent of u>. 

In the strong-coupling regime these modes correspond to the coherent superposition of the cavity field and q s ± 
modes or zero cavity field and q a modes. We can gain further insight about these modes by writing the strong- 
coupling version (a; a "C G) of the equation of motion (|33cp in the sideband-resolved regime 



7 

a a = -^w s a a - ~a a 

7 iG 

a s = -zcj s a s - -a s H -=a 

2 V2 

a = iAa + iV2Ga s — -^a. 



(41) 
(42) 
(43) 



Its eigenfrequencies are given by 



uj a = w s 



lo b± = -A ± 1/4V16G 2 - 7c 2 . 



(44) 
(45) 



Moreover, the eigenvector associated to uja, corresponds (asymptotically) to a vector for which a and a s are zero, 
while the eigenvectors associated to ojb± correspond to vectors for which a a is zero, allowing thus to identify uja with 
the frequency of the dark mode and ujb± to that of the bright modes. 



13 



The coupling G = 7 c /4 corresponds to the onset of the strong-coupling regime where, asymptotically, the cigenmodcs 
correspond to the dark and bright modes. Moreover the limit G>7 C yields 

u B± = -A±G, (46) 

allowing us to identify the symmetric mode q s with the mechanical component of the bright mode, as it can be seen 
from the expression giving the strong-coupling expansion of uj\ (see eqs. (I22p ). 



